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Abstract. Much of our understanding of modern astrophysics rest on the notion 
that the Initial Mass Function (IMF) is universal. Our observations of a sample of HI- 
selected galaxies in the light of Ha and the far-ultraviolet (FUV) challenge this result. 
The flux ratio FhoV/fuv from these star formation tracers shows strong correlations 
with surface-brightness in Ha and the R band: Low Surface Brightness (LSB) galaxies 
have lower Fh^/Zfuv ratios compared to High Surface Brightness (HSB) galaxies as 
well as compared to expectations from equilibrium models of constant star formation 
rate (SFR) using commonly favored IMF parameters. I argue against recent claims 
in the literature that attribute these results to errors in the dust corrections, the micro- 
history of star formation, sample issues or escaping ionizing photons. Instead, the most 
plausible explanation for the correlations is the systematic variations of the upper mass 
limit and/or the slope of the IMF. I present a plausible physical scenario for producing 
the IMF variations, and suggest future research directions. 



1. Introduction 

The assumption that the IMF is constant is commonly adopted in astrophysics. The 
basis for this assumption can be summarized as follows. Color-magnitude diagram 
(CMD) analyses of populous star clusters, where the IMF is e asiest to measure, show 
that the upper end of the IMF has a slope close to the lSalpeterl (11955b value y - -2.35, 
with cluster to cluster variations exp la inable by stoc h astic variations in the number of 



witn cluster to cluster variations expla inable by stoc n astic variations m tne number ot 
stars in each cluster dKroupal 120011) . lLada & Lad"al d2003l hereafter LL03) find that 



almost all star formation occurs in star clusters. Since clusters have a constant IMF and 

all stars form in clusters, then the IMF must be constant. 

■ it 1 1 || I 

Recent studies (including IHoversten & Glazebroo k 2008; Lee et al. 2009) indicate 
problems with this assumption. Our contribution to this field. iMeurer et all (12009. here- 
after M09), has been particularly forthright in its challenge to a constant IMF. Here I 
review results from that study, present a scenario for the origin of the IMF variations 
and respond to critiques of our paper. 



2. The SINGG and SUNGG surveys 



Our results ar e based on two sur veys: the Survey of Ionization in Neutral Gas Galax- 
ies (SINGG: IMeurer et all 120061 hereafter M06) and it's sister Survey of Ultraviolet 
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Figure 1. Our main result: ^hq-Z/fuv correlates with optical surface brightness, 
X. Panels (a), (b) plot the effective surface brightness in Ha and R band respectively. 
Downward and leftward pointing arrows indicate cases with Ha signal to noise ratio 
< 2. The diagonal arrow in the lower right corner of each panel represents the 
average effect of dust absorption, which has been removed. The broken horizontal 
line shows the Fh»//fuv value expected for a constant SFR population having a 
Salpeter IMF extending up to 100 Ai . The curves represent models for a "burst" 
(SFR increase - loops to the right) or "gasp" (SFR decrease - loops to the left) SFH. 
The SFH in each case has a constant base level with the burst or gasp having a 
Gaussian shape with fixed FWHM duration of 10 Myr. Two cases for the ratio of 
maximum to minimum SFR for each SFH are shown 2:1 and 100:1. Additional 
models are shown in M09. 



emission in Neutral Gas Galaxies (SUNGG: IWong 20071). The parent sample of both 
of the se is the HI Parkes All Sky Survey (HIPASS lMever et alj|2004 iKoribalski et all 
120041 ). Nearly 500 galaxies were selected from HIPASS for SINGG so as to evenly 
populate the HI mass function with preferentially the nearest galaxies selected. Over 
300 of these were imaged in Ha and the R band. Nearly 150 galaxies from SINGG 
were selected to be observed by GALEX in the FUV and NUV bands for SUNGG. 104 
galaxies with complete optical and UV datasets were used in M09. 

Our simple HI selection method means we do not suffer any of the common optical 
biases. This is illustrated in M06 where we show that all HI targets in our first data re- 
lease correspond to Ha emitting galaxies (in some cases multiple galaxies ), and that all 



known morphologies of star-forming galaxies are represented. Similarly lHanish et al. 



(2006) demonstrate how our selection allows us to recover the cosmic SFR density of 
the local universe. These results demonstrate that our sample is not biased against any 
particular type of galaxy as long as it has HI. Hence, we can rule out selection effects 
as causing our results; this is an essential component of our claim that we found true 
IMF variations (see 93.2l below). 

The ability to probe the IMF with two star-formation tracers was a major design 
feature of our surveys. Ha emission traces the distribution of the ionizing O stars, 
with initial masses down to At* ~ 2OAl , while the vacuum UV emission detected 
by GALEX traces both O and B stars (initial masses 3Al Q < At* < 20At Q )- The 
total luminosity of a star forming population is dominated by the vacuum UV and is 
fairly insensitive to the parameters governing the upper end of the IMF. The fraction of 
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ionizing light emitted, on the other hand, does depend strongly on the IMF parameters. 
Thus the ratio of Ha to UV emission is very sensitive to the IMF. 



3. The Results 

Our key result is shown in Fig. Q] - the ratio of Ha line flux to FUV flux-density 
Fua //fuv varies strongly with optical surface brightness. In these plots, each point 
represents a single galaxy. The x-axis gives the face-on effective surface brightness X, 
while Fua/fyuY, plotted on the v-axis, corresponds to the ratio integrated over an aper- 
ture large enough to encompass the entire galaxy. By using integrated quantities we 
avoid concerns about small-scale structure and evolution; it doesn't matter that ionizing 
photons may escape HII regions, nor that clusters dissolve, as long as the stars remain 
in the galaxy and the ionization happens in the galaxy. We find that Fh^/Zfuv also 
correlates with other quantities including morphology, luminosity and dynamical mass. 
The lumin osity correlation was independently found in the 11 HUGS sample as dis- 
cusse d bv dLee et al.ll2 009). Since X correlates with luminosity dFerguson & Sandag " 



19881 : iKauffmann et al.ll2003h . it is not surprising that FhoV/fuv correlates with both 



quantities. We find that the correlations with X are stronger and hence will concen- 
trate on them (especially the correlation). The LSB galaxies are most noteworthy 
in these panels. These are typically dwarf irregular (dl) galaxies. Their low FhqV/fuv 
ratios suggest a truncated IMF or steep y. Next, I examine other factors that can affect 
Fua/fpUY- By showing that these can not account for our results, we are left with IMF 
variations as the most likely explanation for the observed correlations. 

3.1. Dust Attenuation 

A key concern with any UV or optical study is dust attenuation. Our corrections for dust 
in the UV are based on the IRX-/3 relation for SUNGG galaxies shown in Fig. [2 k- Our 



calibration falls below the relationship found for starburst galaxies by dMeurer et al. 



and is similar to t he rela tionship f ound i n other samples of local normal galaxies 
(e.g. lSeibert et al.ll2005l: iGil de Paz et al. 2007). In the optical, we employ a secondary 



correlation derived by lHelmboldt et all (12004) from integrated B aimer decrements mea- 
sured in the Nearby Field Galaxy Survey dJansen et ail 2000h . Figure |2p (from M06, 



reproduced by permission of the AAS) shows that our optical dust correction can ac- 
count for the relationship between the ratio of Ha and FIR fluxes and the R luminosity 
for the SINGG galaxies detected by IRAS. Fig.[2]demonstrates that our dust corrections 
can account for the attenuation that ends up in warm dust emission. 

The key arguments against dust causing our correlations are that (a) the typical 
correction ends up to be small compared to the scale of the correlations (as shown by 
the reddening vectors in Fig.Q]), and (b) the correlations between S and Fua/fFVW exist 
even without any dust corrections. This was shown in Fig. 6 of M09, where the LSB 
end of the correlation is significantly below the fiducial Fn a / '/fuv level expected for a 
fully populated Salpeter/Kroupa IMF even before dust correction. Using a larger dust 
corr ection will move the se galaxies further from the fiducial ratio, not closer. 

iBoselli et al.l (120091 . hereafter B09) tested dust corrections similar in form to ours. 
This produced stronger correlations of Fh«//fuv with other quantities in their sample 
than when they performed no corrections. They argued that this is not physically plau- 
sible and hence that such "statistical" corrections are unreliable and should not be used. 
While this is a legitimate concern, their problem does not affect our study. We find 



4 



G.R. Meurer 




Figure 2. Tests of our dust absorption and re-emission model. Panel (a) (left) 
shows the ratio of total far infrared emission to that in the FUV plotted against the 
ultraviolet spectral slope, ft, (derived from t he FUV-NUV color) for the SUNGG 
galaxies also detected in the FIR by IRAS dWongl 12007). Our calibration of the 
global dust reddeni ng rela t ion is shown as the solid line. The calibration for starburst 
galaxies by M eurer et al.l (119991) is shown as the dashed line. Panel (b) shows the 
ratio of Ha and FIR fluxes plotted against R band absolute magnitude. Filled symbols 
mark IRAS detections, while triangles indicate IRAS non-detections as 3<x lower 
limits in Fh»//fir- The curves represent the application of a simple dust reprocessing 
models on stellar population models. The solid line is for a Salpeter IMF over a mass 
range of 1 - 100 M@; the dashed line is for the same y but over the mass range of 1 
to 30 Ai@. The thin solid line and dotted line segment show fits to the data: the thick 
solid line shifted vertically and a simple linear fit respectively. Further details can be 
found in M06. 



that the correlations after dust corrections are somewhat worse than the uncorrected 
correlations, as expected since the corrections are noisy and determined from indepen- 
dent quantities. Our derived average du st attenuation v ector points in nearly the same 
direction as the relationship derived by ICalzettil ( 20011. see Fig. 1 of M09), which is 



satisfying considering how coarse the corrections are. B09 go on to remove sources 
that do not have FIR fluxes from their sample. This results in a much weaker correla- 
tion and very few galaxies with low FhqV/fuv- Dwarf gal axies, typically LSB, have 



very weak FIR emission even in deep Spitzer observations (ILee et al.1 12009: Dal e et al. 



2009). By excluding galaxies without FIR emission, B09 are effectively throwing out 



the most interesting part of their sample. 



3.2. Star Formation History (SFH) 



Rapid c hanges in the SFR ca n lead to large changes in FhqV/fuv- We used Starburst99 
models (Leitherer et al.ll 1999h to test whether SFR variations could cause the observed 



correlations. Our models consist of a population forming stars at a constant rate with 
a Salpeter IMF over the mass range 0.1 to 100 M Q which then experiences a sudden 
increase (burst) or decrease (gasp) in the SFR. Further details of the model are given 
in M09. The results of the modeling are shown in Fig. Q] Panel (a) shows that strong 
changes in the SFR can result in large simultaneous variations in both Fn a //fuv and 
Xho- However, panel (b) shows that the effects on are much smaller; neither the 
beginning of a gasp nor the end of a burst can give the correlated decrease in both 
I,r and FhqV/fuv needed to explain low Fn a /j puv galaxies. Since HI selection does 
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not preclude LSB galaxies with high FhoV/fuv ( nor HSB - low FhoV/fuv galaxies), 
we can not hide galaxies from our selection during inopportune phases of the star for- 
mation cycle. If SFH causes the low FhoV/fuv cases, then all the LSB galaxies in 
nearby universe are simultaneously going through a sharp decrease in the SFR. Such 
synchronicity is implau sible and allows us to rule out S FH as the cause of the observed 
correlations. Similarly, iHoversten & Glazebrookl d2008l) find that the spread of Ha EW 
versus optical colors in SDSS galaxies can only be explained by a bursty SFH if there 
is an implausible synchronicity in the timing of the bursts, and thus ruled out a variable 
SFR as causing their results. 

Boselli (at this conference and in B09) argues that variations in the "micro-history" 
of star formation can explain the large range of Fh«//fuv- However, his toy-model only 
considers the ratio Fh<*//fuv and not surface brightness. This does not address the 
correlations which are at the heart of the problem. Instead he suggests the experiment 
of considering a dwarf galaxy with and without the 30 Dor region. The actual model 
used in B09 consists of regularly spaced bursts with no star formation between bursts. 
The inter-burst period they require to produce the Fu a /fy\jy distribution would also 
produce many Ha non-detections which is inconsistent with SINGG results (M06). 
The dwarf galaxies usually associated with bursts are Blue Compact Dwarfs (BCD). 
Meur er et al.1 d 19991) show that the bursts in BCDs amount to at best a factor of a few 
change in the SFR. Thus an appropriate model for the SFH in dwarf galaxies is the 
factor of two burst/gasp model, shown as the thin lines in Fig. \T\ which are barely 
discernible in these plots because the excursions in FHa//FUV versus X are relatively 
small. Hypothetical extreme bursts are required to explain the Fh«//fuv excursions, 
and even then, the correlation with X# is not explained as noted above. 



3.3. Escaping Ionizing Radiation 

Using Ha as a star forming indicator implicitly assumes that all ionizing photons are 
absorbed by the ISM. The observed Fua/fpvv correlations could occur if the escape 
fraction of ionizing photons, / esc , inversely correlates with E. This scenario would 
require / esc > 0.5 to explain the galaxies in the lowest quartile of Fn a //fuv- While this 
is a large effect, no attempt has been made to observe emission with a rest Ao < 912A 
in LSB galaxies, so it is not yet possible to formally rule out this scenario. Escaping 
ionizing pho t ons a re the preferred explanation of the low Fua/fpvv ratios given by 
Hunter et a l. (2010). They also imply that some Ha is missed from the observations 



because it is below the Eh» detection limit. 

The detection limit concern is easily dismissed. Our measurements and upper 
limits are state of the art and accurate to within the error bars shown, even when Ha 
is not detected in some individual pixels, because those pixels are within the large 
elliptical apertures employed in our analysis (M06). A large / esc also seems unlikely 
for two reasons. First, the problem cases are the LSB galaxies, and these typically have 
a higher ISM content (larger Mhi/Lb) and puffier disks than the larger spiral galaxies. 
Hence, they should be able to retain their ionizing photons better. Second, CMDs of 
the nearest galaxi es indicate that LSB galaxies typically are de ficient in O stars (e.fc 



the nearest galaxies indicate mat LS>c galaxies typically are de ficient in U stars (e.g 
Young et al.ll2007l) compared to HSB irregulars and BCDs (e.g. lAnnibali et"aLll2008h 



Since such studies usually have been on a galaxy per paper basis, then the assumption 
of a fully sampled Salpeter IMF can fool any team of astronomers to find that their 
particular galaxy has a gasping SFH. A better approach would be to obtain uniform 
high quality HST data of many nearby galaxies and then compare CMDs as a function 
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of E. If the sample is large enough, galaxy to galaxy variations will average out and 
an effective IMF can be deduced. The ANGST team have obtained such a dataset 
dDalcanton et alj |2009) and are well poised to do such a study. 



4. What Causes Upper End IMF Variations? 
4.1. Cluster Versus Field Star Formation 

As noted in £l] much of the support for a constant IMF rests on the notion that all star 
formation occurs in star clusters. However, Ha and UV images of galaxies show the 
majority of the light is diffuse rather than being in discernible clusters or HII regions. 
With SINGG we find th at on average ~ 40% of the Ha emission is in HII regions 
while the rest is diffuse (lOey et al.l I2007T) . In UV and U band light, the fraction of 
light that is in comp act clumps (presumably clusters) is typically ~ 20% for starbursts 
(Meu rer et alJll995l) and only a few percent for normal galaxies ( Larsen & Richtlerl 
2000) ; the vast maj ority is diffuse. The dominant diffuse population has low Fh<*//fuv 



dHoopes et al1l200lh and a UV spectrum dominate d by B stars, wh ereas the young "su 



per star clusters" typically have 6 star rich spectra dTremonti et al.ll2001 ). The disparity 



between cluster and field is typically ascribed to infant mortality (e.g. iTremonti et al. 



20011 : iPellerinet all l201Ch . The main problem with this scenario is that naiveley one 



would expect infant mortality to be stronger in the dense environment of st arbursts than 

in nor mal and LSB galaxies, whereas the results of Meure r~et al.l d 19951) and lLarsen & Ric htler 



(2000) indicate that clusters are more prevalent and thus last longer in starbursts. 

The interpretation from LL03 that all star formation occurs in clusters is based on 
an expansive redefinition of the term "cluster" to include objects that are unbound at 
birth (when they shed their natal gas). Better terms for these unbound objects are the 
older terms "group" and "association", which now seem to be largely neglected by the 
extragalactic astronomical community. However, it is important to retain such a dis- 
tinction, because it is crucial in the physics of high mass star formation. Simulations 
indicate that high mass stars can form readily in dense bound c lusters where proto 



2003, 



stars can quickly buil d-up to high mass by comp etitive accretion (Bon nell et al 
l2004h . Alternatively, iKrumholz & McKea (j2008) show that high mass stars can form 
efficiently in a top-down fashion when the ISM density is very high, because then frag- 
mentation is inhibited. In either scenario, the physics is different for star formation in 
dense clumps compared to low density or unbound objects. Recent large area infrared 
surveys reveal a significant "distributed mode" of star formation comprising ~ 3 0% of 
star formation in the solar neighborhood (lAllen et al.ll2007l : iMegeath et al.ll2009 ). This 



star formation, combined with the large fraction of unbound objects that LL03 refer to 
as clusters, is likely to make the majority of the diffuse UV emission in galaxies. 

4.2. Our Scenario 

We posit that the FhqV/fuv ratio variations arise from the difference between star for- 
mation in bound clusters and the field. Bound clusters need a dense cold molecular 
ISM to form. The frac tion of the total ISM that goes into this phase is set by the hy- 
drostatic pressure dMcKee & Qstrikerll 1 977l : iBlitz & Rosolowsky||2006l). which also de- 
termines how tig htly bound clusters are when they form (lElmegreen & Efremovl 119971 : 



Elmegreen 2008). Stars dominate the disk plane potential, thus largely setting the pres- 



sure. This results in the correlation with which also measures the disk plane mass 
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density. A high mass density disk results in high pressure ISM, with a larger fraction of 
its mass in cold dense molecular clouds. These form dense bound clusters that are rich 
in massive O stars, and hence a high FhqV/fuv- Conversely, a low mass density disk 
will have lower X, lower pressure, and less formation of bound clusters relative to field 
stars and unbound objects. This results in a low FhoV/fuv- 



4.3. Implications 

A variable IMF has numerous major implications. Star formation surveys that adopt 
a constant IMF are likely to get a biased result depending on the tracer they use and 
the pressu re of the ISM. U V based SFRs should be more accurate than those derived 
from Ha dLee et al.ll2009h . Similarly, a variable IMF throws in to doubt the basis for 



measuring the SFH from CMDs using the tip of the main sequence as a clock. Finally 
a variable IMF provid es an alternative explanation for the mass metallicity relationship 
dTremonti et al.1 12004) . Instead of requiring a galactic wind to remove excess metals 
from dwarf galaxies, a truncated IMF means that they may not have made them in the 
first place. 



5. Conclusions and Future Work 

The correlation between global Fu a /fy\jy and the surface brightness of galaxies pro- 
vides strong evidence that the IMF is varying, and suggests that the underlying cause is 
the hydrostatic pressure of the ISM which regulates the phase balance of the ISM and 
consequently the efficiency of the highest mass stars to form. As discussed through- 
out this meeting (and noted above) a variable IMF has major implications for many 
branches of astrophysics 

While I believe that the preponderance of evidence clearly demonstrates the vari- 
able nature of the IMF, the astronomical community, as a whole, is not yet convinced. 
One common thread in studies showing IMF variations is the use of Ha as a star forming 
indicator. It would be good to confirm these results using other methods, such as CMD 
studies. Deep CMD studies of large samples of nearby galaxies (e.g. ANGST) are also 
critical for determining how far down in mass the the IMF v ariability goes. Such stud- 



ies would also provide useful tests of the IGIMF hypothesis (Weidner & Kroupall2005l : 



Pflamm-Altenburg et alj|2009h which makes predictions on the form of the galaxy wide 



IMF. Cases of high FhoV/fuv, we U above Salpeter expectations, as seen in our work 
and that of Gunawardhana (this meeting), are also interesting. The highest values may 
be pushing beyond what models can achieve by just adjusting the upper mass limit or 
y, and may instead require truncating the lower end of the IMF at values above 3 vVl Q . 
Such a short-wick IMF is similar to what is required at high redshift to explain dif- 
ferences between the cosmic SFR density and the mass assembly history dDavell2008l : 



Wilkins et al.ll2008bUah . Hence high Fua/frnv galaxies may be good analogs to high 



redshift galaxies. 
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